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Abstract 



We develop the Glauber theory description of initial- and final-state in- 



teractions (IFSI) in quasielastic A(p, 2p) scattering. We study the IFSI- 

H \ 

distortion effects both for the inclusive and exclusive conditions. In inclu- 
sive reaction the important new effect is an interaction between the two 
sets of the trajectories which enter the calculation of IFSI-distorted one- 
body density matrix for inclusive (p, 2p) scattering and are connected with 
incoherent elastic rescatterings of the initial and final protons on spectator 
nucleons. We demonstrate that IFSI-distortions of the missing momentum 
distribution are large over the whole range of missing momentum both for 
inclusive and exclusive reactions and affect in a crucial way the interpreta- 
tion of the BNL data on (p, 2p) scattering. Our numerical results show that 
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in the region of missing momentum p m < 100 — 150 MeV/c the incoherent 
IFSI increase nuclear transparency by 5-10%. The incoherent IFSI become 
dominant at p m > 200 MeV/c. 
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I. INTRODUCTION 



The strength of the initial- and final-state interactions (IFSI) in (p, 2p) scattering is 
usually characterized by the nuclear transparency, T4, defined as a ratio of the exper- 
imentally measured cross section to the theoretical cross section calculated neglecting 
IFSI in the plane wave impulse approximation (PWIA). It is expected that, due to the 
color transparency (CT) phenomenon IFSI effects will vanish and the nuclear trans- 
parency will tend to unity in (p, 2p) reaction in the limit of s — ► 00 and \t\/s ~ 1. /,From 
the point of view of the Glauber-Gribov coupled-channel multiple scattering theory |^rj] 
the vanishing of IFSI corresponds to a cancellation of the rescattering amplitudes with 
elastic (diagonal) and excited (off-diagonal) intermediate states of the initial and final 
protons participating in hard pp scattering. Naive theoretical considerations [[[],[§ suggest 
a monotonic rise of T4 with s in the case of dominance of the point-like perturbative 
mechanism of hard pp scattering ||. However, in the BNL experiment || on large-angle 
(p,2p) scattering at the beam momenta 6-12 GeV/c near 9 c , m , = 90° (here d c .m. is the 
scattering angle in the pp center of mass frame) a decrease of T4 was observed at beam 
momenta > 10 GeV/c. There were suggestions [0-0] that the irregular behavior of Ta 
is due to an interplay of CT effects for hard point-like and non-point-like, resonance or 
Landshoff , mechanisms of large-angle pp scattering. None the less, from our point of view, 
a satisfactory explanation was not found, and up to now the theoretical situation is far 
from being clear. 

In previous works on CT effects in (p, 2p) reaction the IFSI-absorption effects were 
treated within the optical potential approach. This approximation corresponds to taking 
into account only the coherent IFSI. In this case the calculated cross section is related to 
the exclusive (p, 2p) reaction, when the final states of the residual nucleus are exhausted by 
the one-hole excitations of the target nucleus. The allowance for both the coherent and in- 
coherent IFSI corresponds to the inclusive reaction, when all the final states of the residual 
nucleus are involved. The recent Glauber analysis indicates that in the case of (e, e'p) 



scattering the incoherent rescatterings become dominant at high missing momenta (> 250 



MeV/c). Evidently, in (p,2p) scattering, due to the increase of the number of the fast 
protons propagating through the nuclear medium as compared with (e, e'p) reaction, the 
relative effect of the incoherent rescatterings will be enhanced. The theoretical study of the 
inclusive reaction would be of great importance because the data of the BNL experiment 
H correspond namely to the inclusive conditions. The analysis of the missing momentum 
dependence of the nuclear transparency in (p, 2p) reaction within the coupled-channel 
formalism including the incoherent rescatterings invites complications. First, evaluation 
of the contribution of the off-diagonal incoherent rescatterings requires information on 
the off-diagonal resonance-nucleon amplitudes at arbitrary momentum transfer |TI|]I^| . 
Second, inclusion of the incoherent rescatterings make the coupled-channel analysis com- 
plicated from the point of view of the numerical computations. In this situation it is 
reasonable to start the study of IFSI effects in hard (p, 2p) reaction with inclusion of the 
incoherent rescatterings within the one-channel Glauber model. Evidently, only after a 
comparison of the experimental data with the predictions of the Glauber model one can 
understand whether and to which extent the off-diagonal rescatterings or other effects 
are really important. The Glauber analysis [10] of the missing momentum distribution 
in (e, e'p) scattering shows that there is a region of the relatively small missing momenta 
(p m < 150 MeV/c) where the incoherent rescatterings can be neglected. This fact allows 
one to greatly simplify evaluation of CT effects in this region of the missing momentum 



12]. From the point of view of further investigations of CT effects in (p, 2p) scattering it 



is of great importance to clarify whether this is the case in this reaction as well. For the 
above reasons the Glauber analysis of (p, 2p) scattering is highly desirable. In the current 
literature only in ref. Jl3| the Glauber formalism was applied for evaluation of the nuclear 
transparency in inclusive (p, 2p) reaction. However, unjustified approximations made in 
ref. jn| led to a loss of the IFSI-distortion effects (for the criticism of the approach [|T^ 
see ref. [[LTJ and the discussion in section 2 of the present paper). 

In the present paper we evaluate the nuclear transparency in (p, 2p) scattering within 
the Glauber model in the region of p m < 300 MeV/c. In our analysis we neglect the short 



4 



range correlations (SRC) in the target nucleus and describe the nucleus wave function 
within independent particle shell model. In the case of the single particle momentum 



distribution (SPMD) the effects of SRC |14| are still marginal in this region of momenta. 
The analysis of (e,e'p) [p~5|Jl6|1 shows that at p m < 300 MeV/c SRC practically do not 
affect the missing momentum dependence of the nuclear transparency as well. In (p, 2p) 
reaction the distortion effects are enhanced as compared to the case of (e, e'p) scattering. 
It improves the credibility of the independent particle shell model for analysis of (p, 2p) 
scattering as compared to the case of (e, e'p) reaction. 

In our analysis, as in previous works on hard (p, 2p) reaction, we assume the factor- 
ization of hard pp scattering and soft IFSI. We are fully aware that due to the strong 
energy dependence of the cross section of hard pp scattering (oc s -10 ) this approximation 
may be questionable. The qualitative estimates show that the off-shell effects and the 
nonzero energy momentum transfer in soft IFSI in the kinematical region of our interest 
can increase the nuclear transparency by 20-50%. Unfortunately a rigorous evaluation 
of these effects, requiring the relativistic many body approach to the (p, 2p) reaction, is 
not possible at present. However, as we will see the IFSI-distortion effects are typically 
much stronger than the expected magnitude of effects connected with the off-shellness 
of the bound proton and the energy transfer in soft IFSI. Therefore, we believe that the 
factorized approximation for hard pp scattering and soft IFSI is a good starting point for 
evaluation of the IFSI-distortion effects in hard (p, 2p) reaction. 

The paper is organized as follows. In section 2 we set out the Glauber formalism for 
hard (p, 2p) reaction. The numerical results are presented in section 3. In this section 
we also compare the predictions of the Glauber model with the data on the nuclear 
transparency and the missing momentum distribution obtained in the BNL experiment 
PJT7[]. The summary and conclusions are presented in section 4. 
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II. IFSI IN (p, 2p) SCATTERING WITHIN THE GLAUBER FORMALISM 



We begin with the kinematics of hard (p, 2p) reaction. We denote the four momenta 
of the initial and two final protons participating in hard pp scattering as (Ei,pi) and 
(E 3 ,p 3 ), (E^,Pi) respectively. The trajectories of the initial and final protons will be also 
labeled by indexes 1 and 3,4. We use the coordinate system with z- axis along p\ and 
the x-axis in the reaction plane. At high energy for 9 c _ m = 90° the angle between the 
three momenta of the final protons and z-axis in the laboratory frame, Bi,f, , becomes small 

f. ~ [2m p (Ei — m p )] l l 2 / E 1 , here m p is the proton mass). Below we will make use of 
this fact to simplify the numerical calculations. As was stated in section 1, we assume 
the factorization of hard pp scattering and soft IFSI of the fast protons with spectator 
nucleons. Then, the differential cross section of (p, 2p) scattering can be written through 
the distorted spectral function, S(E m ,p m , Ei), in the form 

— (E u t, E m , p m ) = (s, t)S(E m , p m , E x ) , (1) 
ataEmdp^ at 

where da vv j dt is the differential cross section of hard pp scattering, the missing momentum 
and missing energy are defined as p m — p 3 + p^ — p\ , E m — E 3 + E$ — Ei , s is the center 
of mass energy squared of the pp system. To leading order in the missing momentum s is 
given by 

snsJl-ZaA, (2) 

V m v J 

where so = 2m v {E\ + m p ). Notice that keeping in Eq. @ the second order terms in p m 
does not make much sense because the energy momentum transfer in the soft rescatterings 
of the fast initial and final protons in the nuclear medium, which is neglected in the 
factorized approximation ([!]), also gives the effect of the second order in p m . As was 
mentioned in section 1, an accurate treatment of such effects, and of the off-shell effects, 
requires making use of the relativistic many body approach, which goes beyond the scope 
of our exploratory study. 

The distorted spectral function, which, under the factorized approximation ([!]), accu- 
mulates all the IFSI effects, can be written as 



S(E m ,p m ) = \M f (p m )\ 2 5(E m + E A -iipm) ~ m A ) , (3) 
/ 

where Mf(p m ) is the reduced matrix element of the exclusive process p + A4 — > p + p + 
(A — l)f. In Eq. (|3|) and hereafter for the sake of brevity the variable E\ is suppressed. 
In the present paper we will study IFSI-distortion effects at the level of the missing 
momentum distribution, w(p m ), which reads 

w(p m ) = jT^y J dE m S(E m ,p m ) . (4) 
The substitution of (J3|) into yields 

™(Pm) = ^El M /(^)i 2 • (5) 

In our analysis we confine ourselves to a large mass number of the target nucleus 
A> 1. Then, neglecting the center of mass correlations we can write Mf(p m ) as 

M,(M = jfnM;{r 2 , ...,r A )Mn, ...,f A )S(n, ...f^exp^l) . (6) 

Here and ^/ are wave functions of the target and residual nucleus, respectively. The 
nucleon " 1" is chosen to be the struck proton. The use of the same indexes for the spatial 
coordinates of the nucleons in the nucleus wave functions and for labeling the fast protons 
participating in hard pp scattering must not lead to a confusion because we will not use 
the spatial coordinates of the fast initial and final protons. For the sake of brevity, in 
Eq. (^|) and hereafter the spin and isospin variables are suppressed. The factor S(fi, ta) 
in Eq. (||) takes into account the soft IFSI-distortion effects. In the Glauber model it can 
be written in the form 

S(fi, ...,r A ) = Si(fi, ...,r 4 )5 , 3 (fi, ...,r 4 )5 , 4 (fi, ...,f A ) (7) 

where 61,3,4 (r*i, f A ) are the absorptive factors for the initial and final protons which are 
given by 

A 

S i (r 1 ,...,r A ) = l[[l-'yi{r 1 ,f j )} , (8) 



with 



liifiiTj) = 9{ni{rj - r L ))r i (6 i (fi,f j )) . (9) 

Here Hi are the unit vectors defined as n\ = —pi/\Px\ an d ^3,4 = P3a/\Psa\^ bi(ri,fj) = 
[(fj — r*i) 2 — (fii(fj — ri)) 2 ] is the transverse distance between the spectator nucleons 
"j" and the trajectory of the fast (initial or final) proton "i", is the familiar profile 
function of the elastic proton- nucleon scattering (the label "i" reflects the fact that the 
profile function must be calculated at the energy Ei). Eqs. (0), (@) are written under the 
usual assumption that the spectator coordinates can be considered as frozen during prop- 
agation of the fast protons through the nuclear medium. Also, we neglect the interaction 
radius of 90° hard pp scattering, which is expected to be ~ 1/ y/s. 

In our calculations we use for T(b) the standard high-energy parameterization 



r /,v (J tot (pN)(l -ia pN ) 
T{b) = 6XP 



b 



2 



2B pN 



(10) 



Here a P N is the ratio of the real to imaginary part of the forward elastic pN amplitude, 
BpN is the diffractive slope describing the t dependence of the elastic proton-nucleon cross 
section 

da ei (pN) da el (pN) 



dt dt 



exp(-B pN \t\) . (11) 

t=o 

After substituting expression (JD into Eq. (|5]) and making use of the closure relation 

E -,r A W f (r 2 , ...,r A ) = f[ 5{r 3 - r]) (12) 

/ i=2 

to sum over all the final states of the residual nucleus, we can represent w{p m ) in the form 

1 r i i i 

w{p m ) = -j—^ J d 3 fid 3 f[p D (f 1 , r[) exp[ip m (fi - f\)] , (13) 

where 

f A 

PDin,^) = / J| ^^(ri, r 2 , ...,f A )^*(K^2, ■■■,rA)S(fi,r 2 , ...,r A )S*(r[,r 2 , ...,r A ) . 

J j=2 

(14) 



The function p£ ) (f 1 ,f[) may be viewed as an IFSI- modified one-body proton density ma- 
trix. In the PWIA, when the IFSI factors in the right-hand side of Eq. fll4|) equal unity, 
it reduces to the formula for the usual one-body proton density matrix p(fi,f[), and 
Eq. ( |T3D reduces to the expression for SPMD 

1 f l IT ll 

n F(Pm) = ttttq / dfidf[p(fi,f[)exp ip m (ri-f\) ■ ( 15 ) 



(27TJ 

As was stated in section 1, we will describe the target nucleus in the independent 
particle shell model. After neglecting the SRC the A-body semidiagonal density matrix 
^/j(r*i, r 2 , f^)\l/*(f{, r 2 , f^) still contains the Fermi correlations. To carry out the 
integration over the coordinates of the spectator nucleons in Eq. flip) we neglect the 
Fermi correlations and replace the A-body semidiagonal density matrix by the factorized 
form 

A 

^i(ri,f 2 , ...,f^)**(fi,f 2 , ...,r A ) ^ p(r 1 ,r l )Y[p(f i ) . (16) 

i=2 

Here 



z 

is the shell model one-body proton density matrix and <p n are the shell model wave func- 
tions, p A {^) is the normalized to unity nucleon nuclear density. The errors connected 
with ignoring the Fermi correlations must be small because the ratio between the Fermi 
correlation length lp ~ 3/kp and the interaction length corresponding to the interac- 
tion of the fast initial and final protons with the Fermi correlated spectator nucleons 
hnt ~ 4(<7 tot (piV)(72A)) _1 (here (n^) is the average nucleon nuclear density) is a small 
quantity (~ 0.25). Recall, that a high accuracy of the factored approximation for the 
many-body nuclear density in the calculation of the Glauber model attenuation factor for 
the small angle hadron-nucleus scattering is well known for a long time (for an extensive 
review on hA scattering see [|l8|j ). 

After making use of the replacement (|16D in Eq. (|l4l) the missing momentum distri- 
bution (0) can be written as follows 
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(2tt) j 



<pT X <Pi\p(i f \, ri)$(ri, r[) exp[ip m (ri - rj)] , 



where the IFSI factor $(r 1 ,f^) is given by 



f A 

®(ri,r[) = / JJ p A (f j )d 3 f j S(f 1 ,f 2 , ...,f A )S*(r^,r 2 , ■ ■■,r A ) 

J j=2 



d 3 rp A {r)P(r 1 , r 1; r ) 



A-l 



with 



(17) 



18) 



P(ri,f[,f) = [1 - 7i(Fi,f)][l - 7 3 (f 1 ,f)][l -7 4 (fi,f)] 
x[l -7* (fi,f)] [1-7* (f{,f)] [1-7* (f^,f)] . 



(19) 



The right-hand side of Eq. ( Jl9| ) contains the terms up to sixth order in the profile 
functions. To simplify the calculations we neglect the terms which contain the products 
7i(fi, r) 73,4(77, r) and 7*(r l3 r )7 3 4(77, r ). Such terms correspond to simultaneous inter- 
actions of the spectator nucleon with the initial and final protons. Due to the above 
mentioned smallness of the scattering angle for hard pp reaction in the laboratory frame 
(at pi ~ 10 GeV/c 81, /. ~ 25°) they are only important in a narrow vicinity of the spec- 
tator position r with the longitudinal extension considerably smaller than the interaction 
length of fast protons in the nuclear medium. For this reason these terms can be safely 
neglected in calculating the IFSI factor (|18|). It is also worth noting that a rigorous treat- 
ment of such effects requires taking into account the quark content of the proton and can 
not be performed within the Glauber model. 

After neglecting the simultaneous interactions of the spectators with the initial and 
final protons the IFSI factor $(ri,r[) can be written in the form 



$(fi,f£ 



+ e <mrt 

i,j=l,3,4 



1 



1 



A 



£ G«(*)+ £ G\r[) + G M (n) + G 3 \f[] 

i=l,3,4 i=l,3,4 

A-l 

E Gl 4 (n,K)+ E Gr(n,K) + GfM,K) 

i=l,3,4 i=l,3,4 



(20) 



where 
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Gf::i:^K) = (-i) n+m j ^(f)^^,^)...^^,^* (fin-^J^) 

xT il (b il (r 1 ,r))..TabiM,n)rUK(^ (21) 

(here n^(f) = ApA_(r ) is the nuclear density). 

The IFSI factor (|20"D corresponds to the inclusive (p, 2p) reaction, when all the final 
states of the residual nucleus are allowed. In a similar way, starting from the matrix 
element (§) and taking into account in the sum over the final states of the residual nucleus 
in Eq. ([|]) only the one-hole hole excitations of the target nucleus, one can obtain the 
coherent IFSI factor for the exclusive reaction 

® coh (n, K) = S coh (n)S* coh (K) , (22) 



where 




Giin) + Guin) 

i=l,3,4 



A-1 



(23) 



The factorized form of $ co /i(^i, r-J allows one to write the missing momentum distribution 
for exclusive reaction, which we will refer to as w exc (p m ), as a sum of the IFSI-distorted 
distributions for the one-hole excitations 

2 



w exc (p m ) = I d z n(j) n (ri) ey^iiprnr^S^n) 



(24) 



The terms Gi in (|23|) describe the usual attenuation of the initial and final protons in 
the nuclear medium, while the term G 34 is related to the shadowing effect in the system 
of the final protons. It is a connected with the rescatterings of the protons "3" and "4" 
on the same spectator nucleon. The transverse separation of the trajectories "3" and 
"4" is ~ 28i j (z — Z\) (here z is the longitudinal coordinate of the spectator nucleon). 
Hence, the simultaneous interaction of the spectator nucleon with both fast final protons 



is possible up to the longitudinal distance ~ Rint/^i.f. (here Ri nt = y2B(pN) 0.8 fm is 
the interaction radius for soft piV-scattering) from the hard collision vertex. At incident 
beam momentum ~ 10 GeV/c this size becomes as large as the absorption length for the 
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final protons in the nuclear medium. None the less, as we will see in the energy region of 
the BNL experiment || the shadowing correction to the nuclear transparency turns out 
to be relatively small. 

Let us turn to the whole IFSI factor (^0|). The difference between $(fi,f^) and 
^cohifu f[) is connected with the soft incoherent rescatterings of the fast initial and final 
protons in the nuclear medium. The emergence of the interference terms, like G\ with 
i 7^ j, and of the shadowing terms G\ 4 , Gf 4 and G|| in Eq. ( PU[ ) shows that the process of 
excitation of the residual nucleus can not be considered as a simple superposition of the 
effects produced by the initial and final protons. It makes impossible a probabilistic inter- 
pretation of the incoherent IFSI. The incoherent rescatterings make the IFSI factor (|20|) a 
nonfactorized function of f*i and f[. To leading order in 1/A the nonfactorized form of the 
IFSI factor for inclusive reaction is connected with an interaction between the two sets of 
the trajectories (originating from f\ and r[) generated by the overlapping of the functions 
7i(fi, r ) and 7j(r 1; r ) in the integral (|2~lf) for G- functions with i n ,j m > 1- Due to the fast 
decrease of the profile function at b > Rint) the integral (|2l|) , and the interaction between 
the two sets of the trajectories, vanishes unless there is a region where r is close to all 
the trajectories simultaneously which emerge in ([H]) through the functions 7j(fi, r ) and 
7*(r^,r). Simple geometrical consideration shows that the interaction between the two 
sets of the trajectories generated by the incoherent IFSI is for the most part important 
in the region of \y\ — y x \ < Ri n t- In the variables ( ) and (zi — z x ) this interaction 

vanishes more slowly and survive at distances > Ri n t as well. Evidently, as in the case of 
(e, e'p) reaction the short range interaction between the two sets of the trajectories 
at the level of the Glauber absorptive factor will, for the most part, affect the missing mo- 
mentum distribution at p m > 1/Ri n t ~ 200 MeV/c. However, on the contrary to (e, e'p) 
reaction, in (p, 2p) scattering due to the geometry of hard pp scattering, the azimuthal 
symmetry in the (x, y) plane is absent. The effect of the incoherent rescatterings must be 
enhanced in the region of large p m ^ y as compared with the cases of large p m ^ x or p m ,z- 

The numerical calculations can be simplified exploiting the fact that in the integral 
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Ua{t) is a smooth function (in transverse directions) as compared to the profile 
functions. To factor out in (|21|) the transverse integration, we introduce a new coordinate 
system with the z-axis defined as a line which is the center of "gravity" of the trajectories 
ix, i n , ji, j m in the (x,y) plane. Then, in terms of the new longitudinal, £, and 
transverse, r, variables, to leading order in the small parameter Rf nt /R\ ( Ra is the 
nucleus radius) the integral (|2l|) can be written as 



G?C:?£(?ii?i) = J ^nA(O,O,O^(^(r-n))...0(n 5 (f-f 1 ))^(^(f-fJ)^(rr € (f-f 1 )) 

x / d 2 TT il (b h (f h r))..T in (b in (r 1 ,r^^ . (25) 



Here we have written the nuclear density function of the new variables, and the 



vector r must be treated as a vector-function of £ and f. In Eq. (25) we have also taken 
advantage of the smallness of #/./. and made replacements ffjj — > n^, where is the 
unit vector along £-axis. Evidently, similar to neglecting the simultaneous interaction of 
the spectator nucleon with the initial and final fast protons, such a replacement spoils 
the Glauber form of the attenuation factor only in a very narrow region of the spectator 
positions in the vicinity of the hard pp interaction vertex and practically does not affect 



the final numerical results. The integral over the transverse coordinates in Eq. (|25D has 
the Gaussian form and can be calculated analytically. The corresponding formulas are 
somewhat lengthy and we do not present them here. The remaining integration over £ 
was carried out numerically. 

Eqs. (0), (|20D, (|25|) form the basis for evaluation of the missing momentum distri- 
bution (p, 2p) reaction within the Glauber model. Then, the nuclear transparency for a 
certain kinematical domain, D, of the missing momentum can be calculated through the 
formula 

/ d 3 p m w(p 

rn I 

T A (D) = f— -— . (26) 

J d d p m n F {p m ) 

D 

Besides the calculations of the nuclear transparency as a function of p m for the point- 
like domain D T^(p m ) = w(p m )/riF(p m ), in the present paper we calculate the nuclear 
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transparency for the kinematical domains including all the values of the two and three 
components of the missing momentum. The one- dimensional missing momentum distri- 
butions obtained after the integration of w(p m ) over the two components of p m are given 
by 

Wx(p m ,x) = — / dxdx'dydzp(x, y, z, x, y, z)$(x, y, z, x , y, z) exp[ip mtX (x - x)] , (27) 

If, , , , 

Wy(p m ,y) = / dxdydy'dzp(x, y, z, x, y , z)$(x, y, z, x, y , z) exp[ip mty (y - y)) , (28) 

If, , , , 

w z {p m .z) = — / dxdydzdz p(x, y, z, x, y, z)$(x, y, z, x, y, z) exp[ip m . z (z - z')} , (29) 

Z7T J 

The integrated nuclear transparency corresponding to the whole kinematical domain 
of the missing momentum in Eq. (|26| ) is given by 

T A = J dVip A (ri)$(fi,ri)- (30) 
Making use of (|20|) , (|25|) after a simple algebra we can represent (|30D as 

T A = J d*r xPA (r x ) U~\ 



A-l 

8=1,3,4 



(31) 



where 



OC 



i(fi, n,) = / dln A (fi + nil) 
o 

is the partial optical thickness function, and 



<Kn) = G 34 (ri) + G 34 (r*i) + G*(ri, n) + G^ri, ri) + G^n, ri) + G^{n,fi) 

+G 3 M (n, n) + GUn, n) + G$l(?i, n) . (32) 

The terms containing the inelastic pN cross sections in the square brackets in the right- 
hand side of (|3lD describe the usual absorption in propagation of the fast protons in the 
nuclear medium, while 5(fi) yields correction for the shadowing in the final pp-system and 
the interference effects in the incoherent IFSI. Our numerical calculations show that this 
correction is non-negligible. Thus, even in the case of the integrated nuclear transparency, 
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the IFSI do not allow a probabilistic interpretation. However, it is worth noting, that the 
integrated nuclear transparency, as well as the transverse missing momentum distributions 
(E7|), (|28D, are not affected by the interference of amplitudes for rescatterings of the initial 



and final protons. One can see from Eq. (|25|) that the G-functions in the IFSI factor (|20|) 
related to this interference vanish for z% = z'±. 

Eq. (31) yields the nuclear transparency for the inclusive (p, 2p) reaction. In the case 
of the exclusive reaction, after substituting in ( |3T)| ) the coherent IFSI factor (^), one can 
obtain for the integrated nuclear transparency 

A-l 



TT = I d'np^n) 1 - i 



7y H (J m{E 1 ) (1 - ia pN {E l )) t(fi, Hi) - G u {n) 

1 i=l,3,4 



X 1 



1 

A 



£ °lat (Ei) (1 + ia pN (Ei)) t(n, ft,) - GUn) 



i=l,3,4 



A-l 



(33) 



We conclude this section with a short comment on the work [[0| which previously 
considered (p, 2p) scattering within the Glauber model. The authors of [13| also assume 
the factorization of hard pp scattering and soft IFSI, and use the factorized approximation 



(|i6|) for the A-body semidiagonal density matrix. However, then they make use of the 
approximations which can not be justified. First, they neglect in their counterpart of 
our equation (|TJ) the dependence of the IFSI factor on Ar*i = (r*i — r{), and put r*i = 
r 1 = (r*i + r[)/2. Under this approximation for the factorized parametrization of the 
one-body proton density matrix p(r, r) = Pa{\{^ '+ r))W(f— r ) (here W{r — r ) is 
the Fourier transform of the SPMD) they obtained the missing momentum distribution 



which is proportional to the SPMD. Evidently, the approach of ref. [13] misses all the 
distortion effects which, as we shall demonstrate below, are quite strong. Our predictions 
for the integrated nuclear transparency also differ from the results of ref. ||13|| , because 
in the analysis fl3j the shadowing in the final pp-system and interference effects in the 
incoherent IFSI were not taken into account. 
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III. NUMERICAL RESULTS 



In this section we present our numerical results based on the formalism developed 
in the previous section. We performed the calculations for the target nuclei 12 C, 16 0, 
27 Al and 40 Ca. For l2 C and 27 Al we compare the numerical results with the available 
data of the BNL experiment |qjl7|. The calculations were performed for the oscillator 
shell wave functions. We adjusted the oscillator frequency, u osc , for the above set of 
the nuclei to reproduce the experimental values of the root-mean-square radius of the 
charge distribution: (r 2 ) l l 2 = 2.47, 2.73, 3.05 and 3.47 fm for 12 C, 16 0, 27 Al and 40 Ca 
fl9|| , respectively. We obtained the following set of the values of the oscillator radius 
r osc = ( r m p uj osc )~ 1 ^ 2 = 1.59, 1.74, 1.78 and 1.95 fm. The difference between the charge 
distribution and the proton nuclear density connected with the proton charge radius was 
taken into account. The charge density and SPMD in the region of p m < 300 Mev/c, 
calculated with the above set of the oscillator radii are practically indistinguishable from 
the results of more involved Hartree-Fock calculations. The SPMD calculated in the 
harmonic oscillator shell model is also close to the one obtained within a many-body 
approach with realistic nucleon-nucleon potential in ref. J2I| . 

As it was stated in section 2, we use the exponential parameterization of the proton- 
nucleon elastic amplitude. The diffraction slope of the pN scattering was estimated from 
the relation 

g?„,(pAQ(l + aiv) 
N " ItontpN) ' (34) 

In our calculations we define the pN cross sections and a p ^ as mean values of these 
quantities for the pp and pn scattering. We borrowed the experimental data on pp, pn 
cross sections and a pp , a pn from the recent review 

In Figs. 1, 2 and 3 we show the nuclear transparency for w O(p, 2p) and 40 Ca(p,2p) 
reactions at pi^ = 6 and 12 GeV/c (here pi a b = p\ is the incident beam momentum) 
as function of p m for the two transverse and the longitudinal directions of the missing 
momentum, respectively. Besides the results for the whole IFSI factor (p0|) (solid curve) 
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corresponding to the inclusive reaction, we also show the results for the exclusive reaction 
obtained with the coherent IFSI factor fl22|) (dashed curve). To illustrate the role of 
the interference and shadowing effects in the incoherent IFSI we present in Fig. 1-3 the 
results obtained keeping in the square brackets in the right-hand side of Eq. (|20| ) only the 
terms Gi, G l and the diagonal second order terms G\ (dot-dashed curve). In Figs. 1-3 
we also show the nuclear transparency for the exclusive reaction obtained neglecting the 
shadowing term G34 in (^). Figs. 1-3 demonstrate that the IFSI-distortion effects are 
strong both for the inclusive and exclusive reactions. It is seen that the increase of pi a b 
from 6 up to 12 GeV/c practically does not change the nuclear transparency. However, 
as one can see the p m -dependence of the nuclear transparency appears to be sensitive 
to the shell structure of the target nucleus. Figs. 1, 2 show that for the both transverse 
directions of the missing momentum the incoherent IFSI become very important at p m j_ > 
150 — 200 MeV/c. In the case of the parallel kinematics the difference between the nuclear 
transparencies for the inclusive and exclusive reactions is relatively small at \p m ,z\ ^200 
MeV/c. The results presented in Figs. 1-3 indicate that in the case of (p,2p) reaction in 
the region of p m < 100 — 150 MeV/c the incoherent IFSI increase the nuclear transparency 
by 5-10%. The Glauber analysis of (e, e'p) reaction |Tt| yields the relative effect of the 
incoherent rescatterings < 3% in this region of the missing momentum. Thus, we see 
that the transition from (e, e'p) to (p, 2p) reaction enhances the effect of the incoherent 
rescatterings substantially. None the less our analysis indicates that the optical approach 
which neglects the incoherent rescattering is a reasonable starting point for qualitative 
evaluation of the nuclear transparency in (p, 2p) reaction for the relatively small missing 
momenta (p m < 100 — 150 MeV/c). This fact is of great importance for the prospect of the 
further study of CT effects in this reaction. From this point of view it is also important 
that, as one can see from Figs. 1-3, the shadowing correction also does not affect the 
nuclear transparency for exclusive reaction considerably, and can be neglected in the above 
region of the missing momentum. Notice, that Figs. 1-3 demonstrate that neglecting 
the shadowing and interference effects in the incoherent IFSI leads to a considerable 
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underestimate of the contribution of the incoherent rescatterings. 

In Figs. 4-6 we present the nuclear transparency calculated for the kinematical do- 
mains including all the values of the two components of the missing momentum. The 
legend of the curves is the same as in Figs. 1-3. It is seen that the integration over the 
two components of the missing momentum enhances the relative effect of the incoherent 
IFSI as compared with the unintegrated nuclear transparency. The contribution of the 
incoherent rescatterings becomes important even in the case of the p mjZ -dependence of 
T A . 

The results for the integrated nuclear transparency obtained for the same versions 
of the IFSI factors as in Figs. 1-6 are presented in Fig. 7. Fig. 7 further demonstrates 
that neglecting the shadowing and interference terms leads to the underestimate of the 
contribution of the incoherent rescatterings in the inclusive (p, 2p) scattering by the factor 
~ 2. In the case of the exclusive reaction the shadowing in the final pp-system gives rise 
to the increase of Ta by ~ 10%. 

The curves shown in Fig. 6 were obtained according to the formula fl26|) . In ref. || 
the p m>z -dependence of the nuclear transparency were extracted from the cross section 
of (p, 2p) scattering making use of for the one-dimensional SPMD the experimentally 
measured p m . y -distribution (after its normalization to unity). For this reason the data of 
H must be compared with the theoretical nuclear transparency defined as 

rp BNL ( ri \ - Wz (P m > z ^ (oc\ 

A A \Fm,z) norm I \ ' \°°) 

where 

= ^ (36) 

is the normalized to unity one- dimensional IFSI-distorted p m>2/ -distribution. The com- 
parison of the theoretically calculated ratio ( j36|) for the whole IFSI factor ( ^0]) with the 
experimental data of ref. for 27 Al(p,2p) scattering is presented in Fig. 8. As one can 
see, the predictions of the Glauber model are in good agreement with the data || at 
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Plat, = 6 GeV/c. Notice, that the fall down of the theoretically calculated nuclear trans- 
parency with increase of \p m ,z\ is a consequence of the above mentioned enhancement of 
the contribution of the incoherent rescatterings in the region of large \p m ,y\ as compared 
with the case of large |p m ,z|. At pi ao = 10 the experimental values of the transparency in 
the region of p m ^ > are in excess of the theoretical curve by the factor ~ 2. At pi a t = 12 
GeV/c only one experimental point (the bin 200 < p m;Z < 300 MeV/c) overshoots the 
Glauber curve. 

In Fig. 9 we compare the theoretical p mjy -distribution (|36| ) with the BNL experimental 
data for 12 C(p,2p) and 27 Al(j>,2p) reaction at pi a b = 6 and 10 GeV/c. Besides the 
predictions obtained with the whole IFSI factor (solid curve) we show the distribution for 
the exclusive reaction and the one-dimensional SPMD. As one can see the exclusive distri- 
bution and SPMD differ drastically from the experimental distribution. The theoretical 
distribution obtained with the whole IFSI factor is in good agreement with experimental 
data for both nuclei at pi a b = 6 GeV/c. However, at pi a b = 10 GeV/c for 27 Al(p, 2p) reac- 
tion the width of the experimental distribution is in excess of the width of the theoretical 
distribution. In all probability the disagreements of the Glauber model predictions in the 
cases of the p m ^-dependence of the nuclear transparency, Fig. 8, and the p^-distribution, 
Fig. 9, with the data of ref. PJT7|j for 27 Al(p, 2p) reaction at p\ ao = 10 GeV/c are caused 
by the same reason. 



IV. CONCLUSIONS 

The purpose of this work has been a study of IFSI effects in hard (p, 2p) scattering in 
the region of moderate missing momenta \p m \ < 300 MeV/c within the Glauber model. 
To perform such an analysis, we generalized the Glauber theory developed for the small- 
angle hadron-nucleus collisions at high energy, to the case of hard (p, 2p) reaction. We 
studied the missing momentum dependence of IFSI effects both for inclusive and exclusive 
(p, 2p) scattering. The analysis was performed taking into account the shadowing and 
interference effects in IFSI which were not discussed previously. 
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Our numerical results show that the missing momentum distribution in (p, 2p) reaction 
is substantially affected by IFSI as compared to the PWIA case both for the inclusive 
and exclusive conditions. In the inclusive (p, 2p) scattering the incoherent IFSI become 
dominant at p m ± > 150 — 200 MeV/c. Our results show that neglecting the shadowing 
and the interference effects leads to considerable underestimation of the contribution of 
the incoherent rescatterings. In the region of p m < 100 — 150 MeV/c the incoherent IFSI 
increase the missing momentum distribution only by 5-10%. Therefore, at the relatively 
small missing momenta the optical approach is still applicable for a qualitative evaluation 
of IFSI effects. This fact is of great importance for the further investigations of CT 
effects in (p, 2p) reactions because neglecting the incoherent rescatterings simplifies the 
calculations within the coupled-channel approach considerably. Our results also indicate 
that in future experiments it is highly desirable to measure the nuclear transparency 
separately for the regions of small (p m < 150 MeV/c) and high (p m > 200 MeV/c) 
missing momenta. 

For the first time nuclear transparency measured in the BNL experiment || has been 
compared with the theoretically calculated transparency defined according to the prescrip- 
tion of ref. |J. We emphasize that for the strong distortions of W y (p m! y), W z (j> m ,z) an d the 
strong dependence of T^ NL (p m;Z ) on p m<z at fixed beam momentum p lab , it does nor make 
any sense to plot T^ NL vs. p\H = ft a &(l — ^f) an d it is erroneous to interpret the so 
obtained p\H dependence of T^ NL as the true energy dependence of nuclear transparency. 
We also compared the predictions of the Glauber model with the p m ^-distribution ob- 
served in the BNL experiment |17|], such a comparison has also been performed for the 
first time. In both cases a good agreement of the Glauber model predictions with ex- 
periment was found at pi a b = 6 GeV/c. However, the predictions of the Glauber model 
disagree with the data of refs. [0,0 at Pub = 10 GeV/c. Our analysis indicates that the 
discrepancy between the Glauber model results for nuclear transparency and the p m ,y- 
distribution and the data from the BNL experiment at pi a b = 10 GeV/c are likely to be 
connected with the same cause. 
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Finally it is appropriate to comment on the status of predictions of the one-channel 
Glauber model and on the role of the off-diagonal rescatterings in the GeV's energy region. 
There is a widespread opinion that contribution of the off-diagonal rescatterings vanishes 
at the energies much smaller then the CT energy scale E CT ~ RAijnf - m 2 p ) /2 ~ 10 - 20 
GeV (here m* is the mass of the radial excitation of the proton). Our recent analysis [pi 
demonstrates that in the case of (e, e'p) and (p, 2p) reactions this is only the case for the 
coherent rescatterings, while the contribution of the incoherent off-diagonal rescatterings 
becomes small only at the energy of the fast proton (s) E p < E{ 

nc ^ 1-" * p*TTlp*lint ~ 2 3 GeV 

(here li nt ~ {o'tot{pN)(nA)^ 1 is the mean free path of the proton in the nuclear medium, 
T p * is the width of the resonant state). At higher energies (> Eqt) the off-diagonal rescat- 
terings increase the nuclear transparency for the coherent rescatterings and decrease it 
for the incoherent rescatterings. Such an interplay of the coherent and incoherent IFSI 
can lead to an irregular energy dependence of the nuclear transparency in inclusive (p, 2p) 
reaction in the energy region of the BNL experiment. The analysis of ref. [|l^] shows 
that for (e, e'p) reaction the off-diagonal rescatterings may enhance the contribution of 
the incoherent rescatterings by the factor ~ 2 in the energy region of the final proton 
Ei n c < E p < Ect- In the case of (p, 2p) reaction the enhancement factor must be ~ 4. 
Consequently, the off-diagonal incoherent rescatterings may increase the nuclear trans- 
parency in (p, 2p) reaction at pi a b ~ 10 GeV/c by the factor ~ 2. Notice, that correlation 
of the increase of nuclear transparency with the broadening of pm^-distribution observed in 



17fl in "Al(p, 2p) reaction at pi a b = 10 GeV/c gives evidence in favor of the off-diagonal 



incoherent rescatterings cLS cL CcLU.SC of the above rise of nuclear transparency. The above 
estimate corresponds to the off-diagonal rescatterings with the color-singlet intermediate 
3q states, and here it is appropriate to comment on the possible contribution from the 
hidden-color intermediate states. Namely, at high energies the scattering angle in the lab- 
oratory frame 9ij_ 1, and the two 3q states produced in hard pp collision can interact 
with the same spectator nucleon. For this reason the new kind of the hidden-color (for 
instance the (3g){ 8 }(3g){8} state) intermediate states of the two final 3q states produced 
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in hard pp collision may come into play. In the case of the production of the (3g){ 8 }(3g){ 8 } 
state the intermediate hidden-color state can undergo the transition to the normal (3q) (3q) 
state after rescattering on the spectator nucleon (3g){ 8 }(3g){ 8 } + N — > (3g){i}(3g){i} + N. 
It is important that in contrary to the discussed in the present paper shadowing effect in 
the final pp system, the transition of the hidden-color (3g){ 8 }(3g){ 8 } state into the normal 
(3g)(3g) state requires only one Pomeron exchange. Of course, the production amplitude 
for the hidden-color states will be suppressed by the Sudakov form factor. However, in 
the GeV's energy region this mechanism may be potentially important due to the en- 
hancement by the factor ~ N% as compared with the production of the normal (3q)(3q) 
states in hard pp interaction. Thus, we see that one can expect a complicated interplay 
of the diagonal and off-diagonal rescatterings in (p, 2p) reaction in GeV's energy region. 
For this reason, we regard the predictions of the Glauber model only as a baseline which 
allows us to understand the gross features of the coherent and incoherent IFSI and can 
help in further investigation of CT effects in (p, 2p) reaction. 
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Figure captions: 

Fig. 1 - The p mi2: -dependence of nuclear transparency for 16 0(p, 2p) and 40 Ca(p, 2p) scat- 
tering at p mt y = p myZ = 0. The solid curve is for the inclusive reaction, the dashed 
curve is for the exclusive reaction, the dot-dashed curve shows the results obtained 
keeping in the IFSI factor (p0|) only the terms Gi, G l and G\, the dotted curve 
corresponds to the IFSI factor ([22]) without the shadowing term G34 in (E3). 



Fig. 2 - The same as Fig. 1, but for the p miy -dependence of nuclear transparency at 

Pm,z Pm,x 0' 

Fig. 3 - The same as Fig. 1, but for the p m z -dependence of nuclear transparency at 

Pm,x Vm,y 0' 

Fig. 4 - The p m:r -dependence of nuclear transparency integrated over p m y and p m , z - The 
legend of the curves is the same as in Figs. 1-3. 

Fig. 5 - The p m , y - dependence of nuclear transparency integrated over p m z and p m ,x- The 
legend of the curves is the same as in Figs. 1-3. 

Fig. 6 - The p mj2 -dependence of nuclear transparency integrated over p m ^ x and p m ,y The 
legend of the curves is the same as in Figs. 1-3. 

Fig. 7 - Integrated nuclear transparency for 12 C(p, 2p), 16 0(p, 2p), 27 Al(p, 2p) and 40 (p, 2p) 
reactions. The legend of the curves is the same as in Figs. 1-3. 

Fig. 8 - Comparison of the Glauber model predictions for p mi2; -dependence of nuclear 
transparency defined according to Eq. ( pop with the experimental data of ref. for 
27 Al(p,2p) reaction. 

Fig. 9 - The p miy -dependence of the one-dimensional normalized to unity missing momen- 
tum distribution in 12 C(p,2p) and 27 Al(p,2p) reactions. The solid curve is for the 



Glauber model predictions for the whole IFSI factor (|20|), the dashed curve shows 



the predictions of the Glauber model obtained with the coherent IFSI factor ( |22"D 
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the dotted curve shows SPMD, the histogram shows the experimental distribution 



observed in fi% . 
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